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ABSTRACT 

We consider the problem of lensing by binary galaxies idealized as two isothermal spheres. 
Th is is a natural extension o f the problem of lensing by binary point masses first studied 
by lSchneider & We iss ( 1986). In a wide binary, each galaxy possesses individual tangential, 
nearly astroidal, caustics and roundish radial caustics. As the separation of the binary is made 
smaller, the caustics undergo a sequence of metamorphoses. The first metamorphosis occurs 
when the tangential caustics merge to form a single six-cusped caustic, lying interior to the 
radial caustics. At still smaller separations, the six-cusped caustic undergoes the second meta- 
morphosis and splits into a four-cusped caustic and two three-cusped caustics, which shrink 
to zero size (an elliptic umbilic catastrophe) before they enlarge again and move away from 
the origin perpendicular to the binary axis. Finally, a third metamorphosis occurs as the three- 
cusp caustics join the radial caustics, leaving an inner distorted astroid caustic enclosed by 
two outer caustics. The maximum number of images possible is 7. Classifying the multiple 
imaging according to critical isochrones, there are only 8 possibilities: 2 three-image cases, 
3 five-image cases, and 3 seven-image cases. When the isothermal spheres are singular, the 
core images vanish into the central singularity. The number of images may then be 1, 2, 3, 4 
or 5, depending on the source location, and the separation and masses of the pair of lensing 
galaxies. The locations of metamorphoses, and the onset of threefold and fivefold multiple 
imaging, can be worked out analytically in this case. 
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1 INTRODUCTION 

The Cambridge Sloan Survey of Wide Arcs in the Sky (CAS- 
SOWARY, see Belokurov et al. 2007, 2008) has uncovered a num- 
ber of examples of gravitational lensing by luminous red galaxies 
(LRGs). For example, the lens CASSOWARY 2 has two very bright 
and massive LRGs at redshifts of z = 0.426 and 0.432 respectively 
acting as lenses for a blue star-forming galaxy at z = 0.970. Al- 
though the LRGs lie in a loose group, a natural starting point is 
that the lens is a close pair of galaxies. CASSOWARY 5 has two 
early-type galaxies separated by ~ 5" on the sky acting as lenses. 
The three images of the high redshift source visible in the Sloan 
Digital Sky Survey data straddle the two lenses. CASSOWARY'S 
search strategy has proved adept at finding gravitational lenses with 
large separation images (> 3"), which naturally enough are typi- 
cally caused by two or lensing galaxies. There are also a number 
of existing strong lenses for which there ap pears to be more tha n 
one lensing galaxy — such as HE 02 30-213 dWis otzki et al. 1999), 
PMN JO 134-0931 Jwinn et al.ll2003l) and B 1359+ 154 teusin et all 
1200 lh . In fact, in roughly a quarter of all strong lenses, the lensing 
potential is known to be more complex than a single lensing galaxy, 
with usually a pair or group of galaxies being implicated. 

The only binary lens that has so far received detailed scrutiny 
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is the important case of two point masses. [Schneider & Weissl 
( 1986) carried out a detailed study of the caustics, critical curves 
and imaging properties. In particular, they sho wed that the sim- 
pler model of a point ma ss plus external shear JChang & Refsdall 
1 1984 1 An & Evans l2O06h is not always a reliable description of the 
lensing properties of binary point masses. Given the importance 
of microlensing surveys for planets, there has been much s ubse- 
quent theoretical work on the b i nary point mass cas e (see e.g.JWittl 
1 19901 : iMao & Paczynskil Il99ll : Iwitt & Mad 1 19951 : lAsadj l2002h . 
The theory has borne fruit in detailed modelling of bina ry lens 
event s, including the first microlens mass determinations dAn et al.l 
2002), and the discovery of the first Neptune mass extra-solar plan- 
ets teeaulieu et alj|2006t) . 

Here, we provide a theoretical treatment of the problem of 
lensing by two isothermal spheres. The model h as been considered 
numerically in lKochanek & Apostolakisl dl988n . though primarily 
in the context of lensing in two planes at different redshifts. The 
motivation of our study is threefold. First, a binary galaxy model 
has a direct application to the close galaxy pair lenses found by the 
CASSOWARY survey or to groups and clusters acting as lenses. 
Second, the investigation of the properties of highly non-symmetric 
gravitational lenses is seriously incomplete, and so the model pro- 
vides an interesting counterweight to the binary point mass case. 
Third, in the case of binary singular isothermal spheres, many of 
the properties can be found analytically, giving insights that are not 
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so easy to obtain using numerical work. The number of lenses for 
which analytic progress can be made is still very small. 

The paper is organised as follows. In Section [2] we numer- 
ically examine the different possible critical curve topologies of 
two cored isothermal spheres. Different image configurations are 
explored by considering the Fermat surfaces in Section[3] and Sec- 
tion [4] deals with the special case of singular isothermal spheres, 
and derives a number of exact results. 



2 TWO CORED ISOTHERMAL SPHERES 
2.1 Models 

We use dimensionless source pla ne £ = (£, r/) and len s plane co- 
ordinates x = (a;, y) (see e.g., ISchneider et alj[l992h . The two 
isothermal spheres are centred at ±a on the 2>axis. The deflection 
potential is 

ip(x,y) = £1 [r 2 c± + (x + a) 2 + y 2 ] 1/2 

+ E 2 [r 2 C2 + (x - a) 2 + y 2 ] 1/2 , (1) 

where r ci and r C2 are the core radii of the isothermal spheres. If 
the core radii vanish, this reduces to the potential of two singular 
isothermal spheres with Einstein radii E\ and E 2 . We take E\ ^ 
E 2 and a > without loss of generality. 

Let us introduce r\ = (a + x) 2 + y 2 and r 2 — (a — x) 2 + y 2 . 
Then the convergence is 



Ei 2r 2 ci + r 2 E 2 2r 2 2 + r 2 



K ™ 2 (r^+if) 3 /* 1 2 (r 2 2 + r 2 ) 3 / 2 ' 

whilst the shear components are 

E i y 2 -(a + x) 2 E 2 y 2 -(a~x) 2 

711 ' 2 (r^+r 2 ) 3 / 2 2 (r 2 2 + r|)3/2 ' 



(a - x)y 



(a + x)y 



72(X) = (r 2 c 2 + r 2 fli ~ 1 (r\ + r 2 )3/2 ' 
The Jacobian of the lens mapping may be found either from 
det A(x) = [1 - ip xx (x)] [1 - Vro(x)] - V^C*) 2 , 
or from 

det A(x) = [1 - k(x)] 2 - 7 (x) 2 , 
where the magnitude of shear is 



7(x) = \/7i(x) 2 +7 2 (x) 2 



(2) 



(3a) 



(3b) 



(4) 



(5) 



(6) 



For an isolated, singular isothermal sphere, 7 = k and so we obtain 
the familiar result det A = 1 — 2k. The analogue of this in the 
binary, singular case is 



det A = 1 - 2k + AE 1 E 2 



2 2 
a y 



(7) 



and so det A = 1 — 2k on the axis of two singular isothermal 
spheres. 



2.2 Two Identical Galaxies with Varying Separation 

Let us take as properties of a fiducial elliptical galaxy lens a veloc- 
ity dispersion a = 250 kms -1 , a core radius r ci = r C2 = 100 pc, 
and a redshift Z\ = 0.46. Let us take the source redshift as 
z s — 2.15, and use a AC DAI concordance cosmology. As shown 



in IShin & Evanl |2008), this choice of parameters is astrophysi- 
cally reasonable. The corresponding dimensionless Einstein radius 
is Z? « 7.11 (where the length scale is chosen to be lkpc). We 
place two such fiducial isothermal spheres at (—a, 0) and (+a, 0), 
and numerically find critical curves and caustics for varying a. The 
different possible configurations are shown in Fig.Q] 

When a is large (top-most panels of Fig.QJ, the critical curves 
and caustics of the two lenses are disjoint, but the outer tangential 
critical curves are not circular (as they would be for a single such 
lens in isolation) and the corresponding caustics are distorted as- 
troids with four cusps. There are small radial critical curves around 
each lens that map to large roundish caustics. We call this config- 
uration of critical curves, and the corresponding caustic configu- 
ration, 'Type 0'. For a = 10.0, the tangential critical curves have 
merged into a single outer critical curve (which we call the com- 
mon tangential critical curve), and the corresponding caustics have 
merged to form a single six-cusp caustic. We call this arrangement 
'Type 1'. The six-cusp caustic elongates along the y-axis as a de- 
creases further, and two small three-cusp caustics split off, as two 
small critical curves (henceforth 'bean' critical curves) pinch off 
from the outermost critical curve, to leave a single four-cusp astroid 
caustic (a — 5.6 and a = 5.4). We call this 'Type 2'. The physical 
origin of the 'bean' critical curves is explained later (§ 14.4b in the 
context of singular isothermal spheres. At a — 4.0, the three-cusp 
caustics have moved away from the origin and enlarged, whilst the 
'bean' critical curves have moved towards the origin and enlarged. 
Decreasing a still further, the two 'bean' critical curves merge with 
the two radial critical curves, creating two critical curves around 
the origin (bottom-most panels). The inner one maps to a caustic 
with two butterfly cusps, whilst the outer one maps to the outer- 
most caustic. The outermost critical curve still corresponds to an 
astroid caustic around the origin. We call this configuration 'Type 
3b'. 



2.3 Two Galaxies with Different Einstein Radii E\ and E 2 

Although there are five free parameters in (UJ, the geometry of the 
critical curves and caustics are determined by four. This is most 
easily seen from the lens equation 



i = x 



v = y 



£1 



(x + a) 



(r2x+r?) 



2-1V2 



E 2 ■ 



(x- 



Ei 



E 2 



„2\l/2 



(rl 1+ rf) 



2^/2 



+ rlf 2 



(8a) 



(8b) 



If all five parameters are scaled by A, the new lens has the same 
critical curve and caustic geometry, just on a different scale: choos- 
ing new coordinates (x',y') = \~ 1 (x,y), (£',77') = A - 1 77) 
reduces the new lens equation to the old one. 

Now, the most physically interesting regime has r c <a- that 
is, the separation between the two lensing galaxies is much larger 
than their core radii - whilst E\ , E 2 and a are comparable. The 
scaling degeneracy means that it is only the ratios Ei/a and E 2 /a 
that determine the geometry of the critical curves and caustics (up 
to the effect of the core radii). So we fix a (to 10, specifically), al- 
low Ei and E 2 to vary, and find which critical curve geometries 
arise from different parts of the E\-E 2 parameter plane. Since only 
small core radii are physical, we do not explore the r ci and r C2 di- 
mensions of parameter space much. For simplicity, r c \ >c2 are set to 
scale with E\ >2 as (0.1/7.11)^1,2 respectively, so that a doubling 
of both Ei corresponds, physically, to halving the separation 2a be- 



2 



Type 




a = 14.5 


\ ° 




o J- 








-20 -10 





10 20 


Type 1 




a = 10.0 











'. Type 



a = 14.5 






-20 


-10 


10 20 


10 


Type 1 


a = 5.6 




5 








5 















-5 










-10 








-5 




-20 


-10 


10 20 


10 


Type 2 


a = 5.4 




5 




/ ° 




5 















-5 




\ ° 






-10 








-5 




-20 


-10 


10 20 


10 


Type 2 — 


a = 4.0 


5 


5 




/ ^ 























y ^ 






-5 








-10 








-5 



10 20 



-10 


5 10 


Type 1 


a = 10.0 


-10 


5 10 


Type 1 


^^^C^ a = 5.6 

S3) 






-10 


5 10 


Type 2 


^-Z^<7^^ a = 5.4 

(£)) 






-10 


5 10 


Type 2 




" ° = 4 "° 

f ) ) - 





-10 -5 5 10 



Type 3 b 




a = 3.0 



' Type 3b 




a = 3.0 



Figure 1. Critical curves (left panels) and caustics (right panels) for two identical fiducial isothermal spheres (with dimensionless Einstein radius ~ 7.11) 
separated by dimensionless distance 2a. Colours of corresponding critical curves and caustics match. The Type of the critical-curves or caustics is defined in 
the text. 
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tween the lens galaxies. (The factor (0.1/7.11) normalizes the core 
radii to the fiducial isothermal spheres.) 

We find that there are five possible critical curve topologies. 
There are the four seen in § 12.21 Type (disjoint tangential crit- 
ical curves), Type 1 (a common tangential critical curve), Type 2 
(a common tangential critical curve surrounding two 'bean' critical 
curves in addition to the two radial critical curves), and Type 3b 
(the common tangential critical curve, and two other critical curves 
in a distorted annulus around the origin, with the galaxy centres 
lying within the edges of the annulus). In addition, there is a Type 
3a topology: unequal isothermal spheres mean that the two 'bean' 
critical curves merge with one of the radial critical curves before 
the other. These five critical curve topologies are separated by four 
metamorphoses of critical curve geometry, as shown in Figures [2] 
to[5] which we call Metamorphoses 1, 2, 3a and 3b. Metamorphosis 
1 marks the transition from Type to Type 1 critical curves, and so 
on. Recall when viewing these figures that Ei are always dimen- 
sionless Einstein radii (in the same units as the axes) and that the 
isothermal spheres are centred at ±10 on the x-axis. 

The metamorphoses of caustics occur at metamorphoses of 
the critical curves, except for the development of swallowtails in 
one of the caustics without any corresponding change in the topol- 
ogy of the critical curves: see Fig. [6] (For the mathe matical back- 
ground to catastrophe theory, see e.g., lArnoldlll986l whilst for an 
introduction to cr itical-curve and caustic metamorphoses, see e.g. 
iPetters et alj|200lh So we also call caustics Type 0, 1 and so on if 
the critical curves are of that type. Note, however, an astroid caus- 
tic can still pierce a radial caustic, for example, without the critical 
curve topology changing. But caustics do not merge or split unless 
the critical curves do. 

In Type critical curves, det A > outside the tangential 
critical curves, and det A < within them except for the small 
regions within the radial critical curves around the galaxy centres 
at (±a, 0). There is a saddle of det A on the x-axis, and as Ei are 
increased, det A at the saddle decreases. As the saddle decreases 
through zero, the two regions of negative det A link up, forming 
Type 1 critical curves (Metamorphosis 1, Fig[2]l. The corresponding 
metamorphosis in the caustics is the merging of t he two four-cusp 
astroi ds, via a 'beak-to-beak' transition (see e.g.. ISchneider et all 
Il992l Chap. 6), into a six-cusped curve, which we henceforth call 
a 'hexacuspid' . Note also that the hexacuspid can pierce the radial 
caustics. 

As Ei are increased further, the outer, common tangential 
critical curve loses its bottleneck, but eventually develops dimples 
around two small regions that have more positive det A than their 
surroundings (these are regions where the shear is small: see § I4,4| |. 
Meanwhile, the hexacuspid becomes less elongated along the x- 
axis and more elongated in the j/-axis. Increase Ei still more and 
Metamorphosis 2 occurs (Fig.[3j: the dimple in the tangential crit- 
ical curve closes completely around the small regions of positive 
det A, creating two 'bean' critical curves. The two small three- 
cusp caustics split off from the tangential caustic via beak-to-beak 
metamorphoses. 

At still larger Ei, the 'bean' critical curves extend towards the 
radial critical curves whilst the three-cusp caustics move outwards 
towards the radial caustics. When E2 < 20, we have Metamor- 
phosis 3a (Fig.|4] the beans merge with the left-hand radial critical 
curve, forming what we call the 'Pacman' critical curve Q, and the 
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three-cusp caustics tack on to the corresponding radial caustic). As 
Ei reaches < 20, the 'Pacman' critical curve merges with the re- 
maining radial critical curve (Metamorphosis 3b, Fig. [5} to form 
an annular region of positive det A, as the caustics change accord- 
ingly (via hyperbolic umbilic metamorphoses again). 

The five critical curve regimes (Types 0, 1,2, 3a and 3b), into 
which the E1-E2 parameter plane is divided by the four critical 
curve metamorphoses, are shown in Fig. [7] Recall that we have 
fixed a = 10 because the critical curve and caustic geometry is de- 
pendent mainly on the ratios £1,2/(1 and r c 1,2 /a in the astrophys- 
ically important regime of a ^> r c i,2- For small core radii, there is 
a good analytic approximation (which is exact for singular isother- 
mal spheres - see § |4.3t for the lens parameters E\ , Ei , a at which 
Metamorphosis 1 (the merging of the two tangential critical curves) 
occurs. The other metamorphoses are determined by numerically 
finding critical curves and caustics for various Ei. Metamorphoses 
3a and 3b depend very strongly on r c i. C 2, whilst Metamorphoses 1 
and 2 don't, as might be expected. 

Higher order catastrophes such as swa llowtails have been 
noted before in a variety of lens mo dels dKeeton et alj |2000| : 
Evans & Witt 2001; BradaC et al. 2004). Swallowtails can form for 
Types 3a and 3b; regions in the E1-E2 plane where the caustics 
have swallowtails are denoted by 3aS and 3bS. The lower-right of 
Fig- El which is left blank, is - by the symmetry of the lens - just 
the reflection of the upper-left part of the plot in the line E\ = E2 
(shown dotted). 

From Fig. [7] we confirm that Type 0, 1, and 2 curves corre- 
spond to progressively stronger isothermal sphere lenses. The size 
of the core radii affects the size of the radial critical curves and 
hence when the radial and 'bean' critical curves merge via Meta- 
morphoses 3a and 3b. As expected, increasing the core radii mean 
that the mergers take place at smaller Ei, as the radial critical 
curves enlarge (reaching out towards the bean critical curves). The 
reason that smaller core radii enlarge the swallowtail region of pa- 
rameter space is that r c — singular isothermal spheres always 
have three-cusp caustics (see §[4]l. 

It is also worth noting that the 'bean' critical curves and the 
corresponding three-cusp caustics pass through an elliptic umbilic 
catastrophe in parameter space (where they shrink to zero size), in- 
between Metamorphoses 2 and 3a. The condition on Ei, E2 and a 
for this can be found exactly for the case r ci — r C2 — (see § |4.4| l, 
and remains a good approximation for small r c . 



3 IMAGE GEOMETRIES 

There can be 1, 3, 5 or 7 images; the corresponding regions in the 
source plane can be easily deduced in FiguresQ]to[5] Only one im- 
age is produced if the source is outside all the caustics, and there 
are two extra images for each caustic within which the source lies. 
Central images are usually highly demagnified and would not be 
observable, so some 3-image systems would appear as doublets 
(e.g. 3A-1 in Fig. |10t, some 5-image systems as quadruplets (e.g. 
5B-1 in Fig. lilt or even triplets (e.g. 5A-1 in Fig. II It. and 7-image 
systems usually as quintuplets (e.g. 7A-1 in Fig.ll2t. The different 
image geometries may be clas sified according to the topography of 
the Fermat potential (see e.g. JSchneider et al]|l992l) . which is 

€ (x) = i(x-{) 2 -^(x) (9) 

where £ = (^i 7 ?) an d x — ( x ,y)- A helpful visualization of 
eqn. ((9) is provided by the Fermat time delay surface, examples 
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Figure 2. Critical curves (left panels) and caustics (right panels) of Metamorphosis 1 , which takes us from Type (upper panels) to Type 1 (lower panels) 
morphologies. E\ and Ei in the lower panels are greater than or equal to E\ and E2 in the upper panels. As in the case of identical isothermal spheres (Fig.[T}, 
the two tangential critical curves merge, as do the corresponding distorted four-cusp astroid caustics. 
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Figure 3. Critical curves (left panels) and caustics (right panels) of Metamorphosis 2, which takes us from Type 1 (upper panels) to Type 2 (lower panels) 
morphologies. E\ and E2 in the lower panels are greater than or equal to E\ and E2 in the upper panels. Two critical curves pinch off from the outer critical 
curve, corresponding to two three-cusp caustics pinching off from the hexacuspid caustic. 
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Figure 4. Critical curves (left panels) and caustics (right panels) of Metamorphosis '3a'. E\ and E2 in the lower panels are greater than or equal to E\ and 
E2 in the upper panels. (Metamorphoses 3a and 3b are distinct if E\ 7^ £2-) In 3a, the two critical curves corresponding to the three-cusp caustics merge 
with one of the two small 'radial' critical curves to form a 'Pacman' critical curve. The three-cusp caustics merge with the corresponding radial caustic, giving 
it six cusps, although increasing E2 further can reduce that to two cusps via a swallowtail folding (see Fig. [6). 




Figure 5. Critical curves (left panels) and caustics (right panels) of Metamorphosis '3b'. E\ and E2 in the lower panels are greater than or equal to E\ and 
E2 in the upper panels. The 'Pacman' critical curve merges with the remaining radial critical curve; the caustics change accordingly. (Swallowtails - see Fig. [6] 
- are not visible with the curves plotted on this scale.) 
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of which appear below (Figs 1 101 to 1 121). The images occur at sta- 
tionary points of the Fermat surface V</> = 0, which may be local 
maxima, local minima, or saddles. (We exclude cases where y lies 
on a caustic, so there are no images at critical points.) Now, </>(x) 
looks like a concave-up paraboloid at large jx| because tj) — > 
(or a constant), so level curves of 0(x) must be closed. Also, level 
curves cannot cross except at saddle points. A Fermat surface can 
therefore be classified by the topology of the leve l curves that run 
throug h saddles, the critical isochrones, (see e.g. ISchneider et al.l 
( 1992) §5.5). There are four ways to close the critical isochrones 
that run through a saddle, as shown in Fig. [8] which we call Le* 
(the lemniscates) and Li (the limacons) for short. We recall that 
if there are 3, 5 or 7 images, there are 1, 2 or 3 saddle points, re- 
spectively. For a given number of saddle points, it is straightfor- 
ward to count the number of possible critical isochrone topologies: 
it is a matter of counting the number of different ways in which 
limagons and lemniscates can be enclosed within one another con- 
sistently, assuming th at <j>(Sj) ^ 4>(S j) for different saddle points 
St, Sj (seee.g. lSchneider e tal. 1992). There are 2 possible topolo- 



gies for 3-image geometries, 6 for 5-image geometries and 25 for 
7-image ge ometries, but not all of thes e actually arise in a given 
lens model jBlandford & NaravarJI 19861) . 

Although there are many different possible topologies for the 
caustics (Types to 3b, swallowtails, and different overlapping of 
caustics), that divide the source plane up into many different re- 
gions, there are only eight different critical isochrone topologies 
for the double isothermal sphere lens: two three-image cases, three 
five-image cases, and three seven-image cases. These are listed 
in Table [T] and the possible five- and seven-image topologies are 
shown explicitly in Fig. [9] We name them according to the shapes 
of the critical isochrones, starting from the outermost and working 
inwards; the exact placement of the limacons and lemniscates fol- 
lows implicitly. So, for example, Li + Li + Le~ means that an outer 
limacon encloses an inner limacon (and an Li + can only be put in- 
side the minus region of the outer limacon), which in turn contains 
an innermost lemniscate (which must be in the minus region of the 
inner limacon), whilst Li + Le + Le~ means that there is a Le + in 
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Figure 7. The critical curve metamorphoses divide E1-E2 parameter space into regions of different critical curve and caustic geometries. Metamorphoses 1, 
2, 3a and 3b are shown in green, orange, purple and blue, respectively. The swallowtail metamorphosis is shown in red. This plot is for a = 10, which is 
arbitrarily chosen: recall that the critical curves are determined (up to rescaling of axes) by the ratios E±/a, E2/U, r ci /a and r C2 /a. Solid points and lines 
correspond to the fiducial scaling of the core radius r c i,2 = (0.1/7.11)Ei,2 (see i]2.3t , whilst open points are for r c i,2 scaling as 5 X (0.1/7.11)^1,2 and 
plus signs for (1/5) X (0.1/7.11)£a,2. The blue crosses are for Metamorphosis 3b (the metamorphosis most sensitive to core radii) with r c 1,2 fixed instead 
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are very sensitive to them. Increasing r c i 2 suppresses formation of swallowtail cusps and brings metamorphoses 3a and 3b to smaller E±2- The lines of 



metamorphoses 3a and 3b would approach E2 = 2a (= 20 here) and E\ = 2a(= 20) if r c 



0. 



(see 5 143V 



the plus region of the limacon and a Le in the minus region (since 
a Le - cannot be enclosed within a Le + ). 

For the same isochrone topology, image configurations can 
differ qualitatively, as listed in the second column of the Table. So, 
for example, 7B-1 and 7B-2 are two distinct seven-image configu- 
rations corresponding to Li + Le + Le - . Examples of qualitatively 
different image configurations are plotted in Figs [10] to [T2] with 
image positions superimposed on critical curves and Fermat sur- 
faces with critical isochrones. The galaxies are still always centred 
at (±10, 0) in these figures, and the core radii r c i,2 are chosen as 
before. 

Classifying image geometries by critical isochrone topology 
is of more than mathematical interest because it allows the order 
of arrival of images to be deduced (in whole or in part). But, as the 
figures show, systems with the same critical isochrone topology can 
be very different: they can have a different number of core images 
(as in examples 5C-1 and 5C-2 in Fig. II It. or the bright images 
may be positioned very differently (e.g. 7B-1 and 7B-2 in Fig. 112b. 
and, since Fermat surfaces change smoothly with changing lens pa- 
rameters and source position, two systems can have similar image 
positions but have different critical isochrone topologies. Indeed, 



more than one critical isochrone topology is possible for sources in 
the same region of a given critical curve and caustic topology: for 
example, both the Fermat topographies Li + Li + Le - (e.g. 7A-1) 
and Li + Le + Le - (e.g. 7B) can arise from a Type-1 caustic con- 
figuration with the source inside all three caustics (not shown in 
Figures). It is also possible to obtain the same critical isochrone 
topology and similar image positions for lenses with topologically 
distinct critical curves and caustics (e.g. the two 7C-1 panels in 

Fig.rrui. 



4 TWO SINGULAR ISOTHERMAL SPHERES 

When the isothermal spheres are singular (r ci = r C2 — 0), the 
algebra is simpler than the cored (r ci , r C2 > 0) case, and we de- 
rive some analytic results on critical curve and multiple-imaging 
regimes. Specifically, we find the pseudocaustics in § 14.11 which 
allows us (in § |4.2t to derive conditions on the lens parameters 
(Ei , E2, a) for which three- and five-fold imaging becomes pos- 
sible. We find the lens parameters for Metamorphosis 1 (in § 14.3b . 
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Le + Li + 




Le Li 

Figure 8. The four ways to close the critical isochrones of a single saddle point. The left panels show lemniscates, the right panels limacons. Critical isochrone 
topologies for multiple saddles are constructed by combining these basic shapes. 




Li + Le + Le Li + Le + Le + Li + Li + Le 

Figure 9. The five- and seven-image critical isochrone topologies that occur for the double isothermal sphere lens. They are named according to the shape of 
the critical isochrones (see Fig.[8j, starting from the outermost. 
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Figure 10. The different three-image Fermat surface topographies and image configurations. The left-hand panels show critical curves and image positions 
(upper panels of each example), and caustics and the position of the point source (lower panels). Colours of corresponding critical curves and caustics match. 
The right-hand panels show contours of the Fermat potential together with the saddle-points and their critical isochrones (red points and contours) as well as 
the images corresponding to the maximum and the minimum in the Fermat potential (black points). In example 3A-1, one of the three images is a highly- 
demagnified core image, whilst in 3B-1 all three images are of high magnification. The source in 3A-2 is inside a three-cusp caustic but outside both radial 
caustics. The critical isochrone is a limacon, as in 3A-1, but this system would be observed as a triplet: the two images near the 'bean' critical curve are of 
high magnification, and the primary image is also magnified. 

Table 1. The different Fermat surface topographies of a lens consisting of two non-singular isothermal spheres, and the critical curve configurations which can 
give rise to them. 



Number 
of images 



Critical isochrone 
topology 



Name of example 
image configuration 



Critical curve 
and caustic type 



Source positions 



Three images 



Li+ 
Li+ 

Le" 



3A-1 Any type Inside one radial (or 'Pacman') caustic only. 

3A-2 Type 2 Inside a three-cusp caustic only 

3B-1 Type 0, Type 1 Inside astroid or hexacuspid caustic only. 



Five images 



Li+ Li+ 
Li+ Li+ 
Li+ Le" 
Li+ Le+ 
Li+ Le+ 



5A-1 Type 1, 2 Inside both radial caustics only. 

5A-2 Type 2 Inside three-cusp caustic and one radial caustic only. 

5B-1 Any type Inside astroid or hexacuspid and one radial caustic only. 

5C-1 Type 3a(S), 3b(S) Outside the astroid but inside both other caustics. 

5C-2 Type 3aS Inside swallowtail but outside radial caustic. 



Seven images Li + Li + Le 
Li+ Le+ Le" 



Li+ Le+ Le" 
Li+ Le+ Le+ 



7A-1 
7B-1 



7B-2 
7C-1 



Type 1 
Any type 



Type 2 
Type2 

Types 3aS, 3bS 



Inside hexacuspid and both radial caustics. 

Type 1 : inside hexacuspid and both radial caustics. 

Type 2: inside both radial caustics and the astroid or a three-cusp caustic. 

Types 3a(S), 3b(S): inside astroid and both other caustics. 

Inside both radial caustics and a three-cusp caustic. 

Inside both radial caustics and a three-cusp caustic. 

Inside swallowtail and the radial caustic. 
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5A-1 




-30 -20 -10 10 20 30 



Figure 11. The different five-image Fermat surface topographies and image configurations (left- and right-hand panels as for Fig. llOt . The saddle-point images 
coloured red arrive before those coloured green. Examples 5A-1 and 5C-1 would be observed as triplets because of the two highly-demagnified images near 
the galaxy centres at (±10, 0), whilst the others would be observed as quadruplets. 5A-2 and 5C-2 are qualitatively different image configurations that have 
the same critical isochrone topologies as 5A-1 and 5C-1, respectively. In 5A-2, the source is placed inside a three-cusp caustic and one of the radial critical 
curves, so only one of the five images is a highly demagnified core image. In 5C-2, a different choice of r ci and r C2 is made so that a swallowtail cusp in the 
caustic that corresponds to the 'Pacman' critical curve can pierce the radial caustic. A source placed inside the swallowtail but outside the radial caustic gives 
five images, of which only one is highly demagnified. 



the elliptic umbilic catastrophe in Type 2 critical curves (§ !4.4b . and 
for the singular limit of Metamorphoses 3a and 3b (§E3}. 



4.1 Pseudocaustics 

The critical curves and caustics can still be of Type 0, 1 and 2 
(with a mild abuse of terminology). The difference between the 
singular and cored cases is that the radial critical curves around the 
two galaxy centres shrink to zero size as the core radii decrease 
to zero. This means that there are no radial caustics, and that any 
core images move towards the centres of the isothermal spheres 
as r ci , r C2 — > 0, their magnifications diminishing, until they dis- 
appear at the singularities. A source can therefore produce either 
an odd or an even number of images, depending on the number of 



core images it would produce if the core radii were not zero. The 
number of images changes by one (with an image being created 
or destroyed at the centre of a singular i sothermal sphere) when th e 
source crosses a pseudocaustic (see e.g. lEvans & Wilkinsonl 19981) . 
a closed curve in the source plane which is approached by the radial 
caustics as r ci , r C2 — > 0. 

The positions of the pseudocaustics are easily found by con- 
sidering the lens mapping on small circles around the singularities. 
In vector form, the lens equation reads 

£ = x - E 1 ^±?- - £2^-^, (10) 
x + a x - a 

where a = (a, 0). Consider a small circle, say of radius S, in the 
lens plane centred at the second singular isothermal sphere at a. 
Then, as 5 — > 0, then x — *• a and (x+a) — > (2a, 0), whilst E2 (x— 
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Figure 12. The different seven-image Fermat surface topographies and image configurations. The order of arrival of saddle-point images is red, green, blue. 
The two 7C-1 systems are not qualitatively different (although the image near (—10, 0) is more demagnified), but arise from a topologically distinct caustic 
configurations. Except for the 7C-1 from Type 3aS caustics, and possibly the 7C-1 from Type 2 caustics, these would all be observed as quintuplets, as they 
each have two core images. 7B-2 has the same critical isochrone topology as 7B-1, but the image positions are very different and one of its critical isochrones 
(coloured red) is topologically a lemniscate but extremely distorted. 



a) / |x — a| becomes E2 times a unit direction vector. So, the circle 
around the singularity at a maps to a circle in the source plane 
centred at (a — E± , 0) and with radius Ei . The other pseudocaustic, 
corresponding to the singularity at —a is found in the same way to 
be circle of radius E\, centred at (— a + E2, 0). 

Figure [T3l shows the critical curve and caustic configuration 
for a pair of singular isothermal spheres that, if cores were put 
in, would have been of Type 2 and near Metamorphosis 3a. The 
number of images produced by point sources in different regions 
is shown on the plot of the caustics. Since the radial critical curves 
have shrunk to points, Type 3a and 3b critical curves do not arise 
in the case of singular isothermal spheres. Instead, there remain 
two three-cusp caustics corresponding to the two small 'bean' crit- 
ical curves of the Type 2 configuration, although these caustics be- 
come indefinitely elongated along one (both) pseudocaustics for 
E2 ^ 2a (respectively, E\ ^ 2a), as the 'bean' critical curves then 
touch the singularities at (±a, 0) (see § [43}. 



4.2 The Onset of Three- and Five-Fold imaging 

When the singular isothermal spheres are far apart, so that their 
caustics and pseudocaustics do not overlap or join, then a point 
source can produce 1, 2 or 4 images. Each isothermal sphere has 
a circular pseudocaustic and a small astroid caustic within it, the 
size of which shrinks to zero as a — > 00. A source outside the 
pseudocaustics produces 1 image, a source inside one of the pseu- 
docaustics produces 2, and a source inside the astroid produces 4. 

As a is decreased, triple imaging becomes possible when the 
pseudocaustics overlap. The condition on the lens parameters for 
overlapping pseudocaustics follows from their positions. It is 

E 1 +E 2 >a. (11) 

Five-fold imaging becomes possible once part of an astroid 
caustic lies inside both pseudocaustics. The pseudocaustics are 
given in § 14.11 and it is straightforward to find the x-intercepts of 
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Figure 13. The left panel shows critical curves, and the right panel caustics (solid lines), for two singular isothermal spheres at (—10, 0) and (+10, 0) (marked 
by magenta and ochre dots). The pseudocaustics are the dotted circles; when a source crosses a pseudocaustic, a single image is created or destroyed at the 
corresponding singularity. The number of images for sources in different regions of the source plane are shown. If r e , were raised from zero, the magenta 
pseudocaustic would become a true caustic as a small radial critical curve formed around (—10, 0), and if r ci were increased further, the three-cusp caustics 
would join with the magenta caustic as a 'Pacman' critical curve formed (Metamorphosis 3a). Raising r C2 from zero would make the ochre pseudocaustic a 
true caustic. 




Figure 14. Close-up of the piercing of a pseudocaustic by an astroid caustic at the onset of five-fold imaging by a double-SIS lens. Ei increase from left to 
right panels. The astroid caustics are within their respective pseudocaustics unless the pseudocaustics overlap, when the astroid of the less-massive SIS (taken 
w.l.o.g. to be the left-hand one) pierces the pseudocaustic of that SIS, and therefore also overlaps with the region enclosed by both pseudocaustics. 



the tangential critical curve(s), and hence the ^-intercepts of the as- 
troid caustic(s). This allows us to find the conditions on the lens 
parameters for which the pseudocaustics overlap with the astroid. 
Such overlapping is possible not only for Type 1 and 2 caustics, 
but also for Type 0. It turns out that, as a is decreased (or Ei in- 
creased), five-fold imaging becomes possible exactly when triple 
imaging does: when critical curves are of Type 0, the astroid of 
the second isothermal sphere stays wholly within its pseudocaustic 
(note we still take Ei < E2, without loss of generality), whilst the 
astroid of the first isothermal sphere pierces its pseudocaustic if and 
only if the pseudocaustics overlap (see Fig. 1141). We now show this. 



The critical curve of the first isothermal sphere cuts the x-axis 
at two points, and it is the right-hand one - call it xn - that maps 
to the right-hand cusp of the astroid that pierces the pseudocaustic 



as shown in Fig. [14] From equations ((2} and Q, the Jacobian along 



the a>axis for r c 



det A(x,0) = 1 



Si 



is simple: 

E-2 



\x + a \x — a 



(12) 



Now, xti lies between —a and a, and the appropriate root of the 
quadratic dl2t is 

x tl = ~ (E! - E 2 - VB 2 - 4C) , (13) 

where B — E2 — E\ and C = a(Ei — E2) — a 2 . The lens equation 
gives 

£« = x tl -E 1 + E 2 . (14) 

Call the right-hand ^-intercept of the first pseudocaustic £ ps i- We 
know from W.llthat 
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-a + Eh + Ei. 



(15) 



Now, the astroid pierces the pseudocaustic if and only if £u > £ ps i, 
that is, combining the three equations above, if 



2a-3E 1 -E 2 -s/(E 2 - E^ 2 + 4a(a - E x - E 2 ) > 0,(16) 

which, on squaring and cancelling terms, is equivalent to 

(a - 2E 1 ){a - Ek - E 2 ) > a{a -Ex- E 2 ) . (17) 

We are considering critical curves of Type 0, and so - anticipating 
the result j23\ from the next section together with the assumption 
Ei ^ E 2 - 03 reduces to 



a - Ei - E 2 < , 



(18) 



which is exactly the condition il lb for overlapping pseudocaustics. 
(A similar calculation shows that the astroid caustic of the second 
isothermal sphere does not pierce the second pseudocaustic.) 

Calculations for the left- and right-most astroid cusps on the 
£-axis, corresponding to critical points at x < —a and x > a, show 
that five-fold imaging is always possible if the critical curves and 
caustics are of Type 1 or Type 2 (with E\ ^ E 2 ^ 2a). However, 
the astroid may lie entirely outside the first pseudocaustic if E 2 > 
2a. In this case, there is still overlapping if 



E 2 (E 2 - Ei) < a(2E 2 - E x ) 



(19) 



4.3 The First Metamorphosis 

By © and direct differentiation of the deflection potential (TJ, the 
Jacobian of the lens mapping for two isothermal spheres reduces, 
along the :r-axis, to 



det A(a;,0) = R(x)T(x) , 
where roots of 



R{x) = 1 
T(x) = 1 



Eir 2 ci 



E 2 r 2 c 



[r 2 ci +(x + a) 2 ] 3 ? 2 
Ei 



+ (x — a) 
Ei 



213/2 ■ 



[r% + (x + a)2] 1/2 [r% + (x~ a) 2 ] 1 ' 2 : 



(20) 



(21a) 



(21b) 



give radial and tangential critical points, respectively. If T(x) = 
has two roots for x £ (—a, a), then the critical curves are of Type 
(there are two separate tangential critical curves). For singular 
isothermal spheres, R(x) = 1. There are no radial critical curves, 
and the three equations above reduce to J 1 2b . that is, 

E 2 



det A{x,0) = T{x) 



Ei 



For x g (—a, a), T(x) is concave down with a single maximum 
at some x m . T'(x) = is a quadratic in x with one root in this 
interval, and we find 



fE 2 



/El 



(22) 



VE 2 + s/Ei 

If T(xm) > 0, there are two roots; if T(x m ) < 0, there are no 
roots. If T(x,„) — 0, the tangential critical curves are touching and 
this is Metamorphosis 1. Combining the two previous equations, 
T(x m ) reads 

J_ , ^2 

2a~ 



1 



(23) 



Critical curves are of Type if the left-hand side (LHS) is less 
than 1 and Type 1 if LHS > 1. It is clear that d23t remains a 



very good approximation for the non-singular isothermal spheres 

if r c , , r c , <g; a. 



4.4 The Elliptic Umbilic Catastrophe 

The radial and tangential critical curves are well-known from sin- 
gle spherical or elliptical lenses. In such cases, the Jacobian of the 
lens mapping, det A, is positive immediately around the galaxy 
centre, negative between the radial and tangential critical curves, 
and positive outside the tangential critical curve, where det A — > 1 
as |x| — 00. We intuitively expect a double isothermal sphere lens 
to have Type critical curves when the separation 2a is large, and 
for the tangential critical curves to merge to give Type 1 when the 
isothermal spheres are brought closer together. The possibility of 
Type 2 critical curves, with 'bean' critical curves, may not be as 
intuitively obvious. 

The 'bean' critical curves surround small regions of positive 
det A inside the common tangential critical curve, where det A 
is otherwise negative. The physical reason that these regions must 
exist is that there are points where the shear is zer o (c.f. discus- 
sion of the two point-mass lens in ISchneider et af] Jl992h §8.3). 
From the expression (O for det A, if 7(x) = at some x, then 
det A(x) > 0. 

For the double isothermal sphere lens, there are indeed two 
points where 7 = because the shears from the two isothermal 
spheres cancel out. Let us begin by showing that, if the isothermal 
spheres (singular or not) are located at (±a, 0), the points of zero 
shear lie on a circle of radius a centred at the origin. It is useful to 
define the angle of shear u> in the usual manner and write the shear 
as the complex (polar) quantity 

2iu 



71 + 172 = 7 e 



(24) 



(Note that w w w ± jt has no physical effect.) Now, since the 
deflection potential of two isothermal spheres is the sum of their 
two separate deflection potentials, the (complex) shear is also the 
sum of two separate shears: Y = L' 1 ' +T^ 2 \ Therefore, if the shear 
is zero at some point x, we must have 

(2k + 1)tt 



Mi 



k € 



(25) 



For an isothermal sphere, the direction of shear uj^\x) at a point 
x must, by symmetry, be either 



or ^(x) = 0j — 



2 ' 



(26) 



where 9j is the direction of x from the centre of the jth isothermal 
sphere. (It turns out to be the latter.) Therefore, l |25t holds if 



0a - 0i = ±2 



(27) 



that is, if the line segment joining (—a, 0) to x is perpendicular to 
that joining (a, 0) to x. The locus of such points x is the circle with 
the two galaxy centres as its diameter, |r| — a (see Fig.|15b. 

It remains to show that there must exist two points on this 
circle for which 7' 1 ) = For a single isothermal sphere, the 
magnitude of shear is 



v (i) 



'-Ei 



( r 2 +r 2)(3/2) ' 



(28) 



where tj is the distance from its centre. As we are only considering 
r cl ,r C2 <SC a,Ei,E 2 , at points on the circle close enough to 
the first isothermal sphere (the one at (—a, 0)), we have 7' 1 ' > 
7^ 2 ', whilst the reverse inequality holds at points close enough to 
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Figure 15. The point of zero shear lies on the circle |r| = a. The centres of 
the (singular or cored) isothermal spheres are at ±a on the x-axis. 



the second isothermal sphere. There are therefore two points in- 
between (one in the y > semicircle and the other in the y < 0) 
where equality holds. 

Although this result does not allow the conditions for Meta- 
morphosis 2 to be determined analytically, it does allow us to find 
the lens parameters at which the elliptic umbilic catastrophe oc- 
curs: it occurs when the point of zero shear is also a critical point. 
We first need to know the (polar) coordinates (a, 9 Z3 ) of the point 
of zero shear. Recalling that n and ri are the distances between 
a point r and the centres of the isothermal spheres, on the circle 
Irl = a we have 



rl = 2a 2 (l +cosf 



r 2 2 = 2a 2 (l 



(29) 



by the cosine rule (see Fig. 115b - Now, by {28}, the point of zero 
shear for the double lens has 

n = Ei 

and combining this with j29\ yields the result 



(30) 



COS f zs = 



El - El 



(31) 



Ef + Ei 

An elliptic umbilic catastrophe occurs at the point of zero 
shear (where the surrounding critical curves - and the correspond- 
ing three-cusp caustics - shrink to zero size) when 



(32) 



This is proved as follows. From {5]l, at the point of zero shear, 
detA = (l- K ) 2 (33) 
where the convergence k is 

K ( r ,0) = K « W 2) = — + — . (34) 
2ri 2r2 

Now, r\ and V2 at the point of zero shear are given in terms of 
Ex, E 2 and a by l|29j and OT}, and putting OU into (1 - k) 2 = 
yields, after some algebra, < |3 2 b - 



4.5 The Singular Limit of Metamorphoses 3a and b 

As E2 — > 2a or Ei — > 2a from below, the 'bean' critical curves 
(that surround points of zero shear) approach the singularity at 
(— a, 0) or (a,0), respectively, as illustrated in Fig. [T6] This can 



be seen from two observations. 

First, near the singularities at (±a,0), det A < on the x- 
axis by inspection of d 12b - Secondly, along the circle [r| = a, 
det A < near (a, 0) or (— a, 0) if and only if E\ < 2a or 
E2 < 2a, respectively. This follows from inspection of 



detA(a,9) = 



1 



E1E2 + 2a z |sin( 



2a 2 sin 6 
2a (Ex |sin(0/2)| + E 2 |cos(0/2)| 



(35) 



which can be obtained - a little tediously but straightforwardly 
either from 



det A(r,8) = - 



d^/dr di/dO 
drj/dr d-q/dO 



(36) 



(with the derivatives being evaluated directly from the lens equation 
in polar coordinates), or from <[5jl , {29), l |34t and 



7M) 



„(i) 



„(2) 



rx 



£2 
r 2 



(37) 



5 EXTERNAL SHEAR OR INTERNAL ELLIPTICITY 

In this section, we investigate two astrophysically important effects. 
First, Binary galaxies may be embedded in a loose group or cluster 
(as in fact is the case for CASSOWARY 2), in which case tidal 
shear will distort the caustics and critical curves. Second, galaxies 
are characteristically flattened, and so internal ellipticity will also 
play a role. In this section, we introduce simple models in § 15.11 
survey the principal effects in § 15. 21 and give some analytic results 
in SO 



5.1 Models 

5.1.1 External Shear 

The effects of a larger-scale gravitational field, arising from the en- 
vironment in which the lens sits, are often approximated by external 
shear: the quadratic terms of the large-scale field expanded about 
the centre of the lens. With suitable translation of the source-plane 
coordinates, the lens equation including external shear is 



i>x(x,y) 



(38) 



where ip is the deflection potential of the two isothermal spheres 
and 
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7 cos 2(j> 7 sin 2<f> 
7 sin 2<j> —7 cos 2<j> 



(39) 



is the external shear matrix, and 7 and (f> are the magnitude and 
angle of external shear respectively. (Note that <f> — corresponds 
to shear from a perturbing mass at 00 or —00 on the y-axis). The 
Jacobian of the lens mapping becomes 

det A = (1 - V** ~ Pu) (1 - i>vv ~ P22)-(i> xy + P12) 2 -(40) 

By the arguments of §|4] the pseudocaustics are still circles of radii 
Ei and E2, but are now centred at 

(— a + Ei + 7a cos 2<j> , ya sin 2<j>) (41a) 
and 

(a — E\ — 7a cos 2(j> , — 7a sin 20) . (41b) 
5.7.2 Ellipticity 

Galaxies are also characteristically flattened. One commonly used 
flattened generalization of the isother mal sphere is the pseudo- 
isothermal elliptic po t ential (see.e.g . Kassiola & Kovneri 1 19931 ; 
iHunter & Evansl l200ll ; lEva ns & Hun terl 120021) . which coincides 
with the isothermal sphere model when its ellipticity parameter e 
is zero, and otherwise has elliptical isopotentials with axis ratio 
(1 + e)/(l — e). T he models are the projections of the power-law 
galaxies studied bv lEvansl fl994). When the isothermal spheres at 
(±a, 0) are replaced with this model, the deflection potential of the 
lens becomes 



where 



Ei 



E2 



(r 2 +r 2 )V 2 (r 2 +r 2y/2 



ri = x' 2 (l-e)+y' 2 (l + e) 
r\ = x " 2 (l-e)+y" 2 (l + e) 



(42) 



(43) 



and (x',y') and (x",y") are Cartesians centred on (— a,0) and 
(a, 0) inclined at <j>i and (j>2 to the x-axis respectively. 



5.2 Effects on Critical Curves, Caustics and Image 
Multiplicity 

A numerical investigation of the critical curves and caustics illus- 
trates the essential differences between these models and the two 
isothermal spheres. For realistic values of shear or ellipticity, no 
dramatically different critical curve topologies arise, although there 
are minor differences resulting from the breaking of the reflection 
symmetry about the x-axis: Metamorphoses 2, 3a, and 3b can occur 
in two steps, with one 'bean' critical curve being created, or merg- 
ing with a radial critical curve, before the other one - examples are 
presented in the panels of Figs.|17|and|18l 

More interestingly perhaps, the maximum number of images 
can change with the addition of shear or ellipticity. Of course, when 
shear or ellipticity is added to a single, isolated isothermal sphere, 
the maximum possible number of images increases from three to 
five (or from two to four in the singular case) because the tangential 
caustic, which had been a single point, becomes an astroid caustic 
with the breaking of circular symmetry. There is not always such 
a dramatic difference when shear or ellipticity is added to a dou- 
ble isothermal sphere lens, because circular symmetry has already 
been broken. However, the tangential critical curve(s) is (are) dis- 
torted by shear (see, for example, the left panel of Fig. 121 1. Shear 



with <f) ~ 0°, or ellipticity with cj>i ~ 90°, can lead to extra cusps 
forming in the tangential caustic(s), which can then self-overlap. 
When there are two separate tangential caustics, this can increase 
the maximum number of images from five to seven (from four to 
six in the singular case); when there is a common tangential caustic, 
there can be regions of nine- and even eleven-fold imaging (seven- 
or nine-fold for the singular case). 

The left panel of Fig. [19] shows the caustics for two cored 
isothermal spheres in moderately strong external shear aligned 
close to the x-axis. The originally six-cusped caustic develops four 
extra cusps, and the self-overlapping creates regions of high-order 
imaging: Source 1 is imaged nine times, and Source 2 imaged 
eleven times. The right-hand panels of Fig[l9] show the Fermat 
surfaces for the two sources, along with image positions and the 
critical isochrones. Note that the elevenfold imaging configuration 
bears a remarkable similarity to the recently discovered lens CAS- 
SOWARY 2, although of course this is not a fit. 



5.3 Analytic Results 

One simple result remains analytic for a double isothermal sphere 
lens in external shear: the touching of the two critical curves along 
the x-axis. When the isothermal spheres are far apart, the shear 
distorts the tangential caustics into ellipses with minor axes aligned 
with the direction of external shear, as ex pected from the effect 
of shear on one isothermal sphere (see e.g. [Schneider et alj|2006l) . 
The tangential caustics are of course no longer ellipses when the 
isothermal spheres are brought closer together, but they still tend to 
be stretched in the direction perpendicular to the external shear. We 
expect shear with (j> = 0° (90°) to shift Metamorphosis 1 to larger 
(smaller) E\ , E2 compared to the shearless case. This is indeed the 
case. 

There is a complication, however. Since even a small nonzero 
7 at (f> ^ 0°, 90° distorts the critical curves (see Figs. 1201 and!21t 
and destroys the reflection symmetry in the x-axis, it is no longer 
clear that det A(x, y) must have a saddle point on the a>axis as it 
did when there was no external shear. Indeed, the two tangential 
critical curves do not necessarily first merge at the a>axis: for very 
strong shear (7 > 0.45 for a — 10) and small <f>, the tangential 
critical curves may merge at two points above and below the x- 
axis first, giving rise to new critical curve geometries, examples 
of which are shown in the panels of Fig. [20] These new critical 
curve metamorphoses don't occur if there is internal ellipticity but 
no shear. 

However, as we now show, in the case of singular isothermal 
spheres, there is still always a saddle point in det A(x,y) on the 
x-axis, which, for reasonable values of 7, is the point at which the 
tangential critical curves first touch (when det A = at the saddle 
point). First, we obtain V det A from d40t and direct differentia- 
tion. Restricted to the x-axis, we have 



0_ 

dx 



det A = 



Ei (x + a) E2 (x — a) 



(1-7 cos 20) , (44a) 



| detA =f^±# + ^#) >-„„2 () . ,44b. 



x + a ' 



\x — a\ 



so it is clear that there is one stationary point on the x-axis for x £ 
(—a, a) because det A| v =q and det A| a =o have a common 
root. This stationary point is a maximum of 

£1 



detA(x,0) = l + 7cos20 



E2 



x + a\ x — a\ ^ 
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x (1 — 7 cos 20) — 7 sin 2<f> . 



(45) 



which for 7 = reduces to d!2t . By the argument of § 14.31 the 
tangential critical curves touch each other on the x-axis when 



(1 — 7 cos 2cf>) (1 + 7 cos 2(f> — 



2a 



— 7 2 sin 2 2(j> 



(46) 



The condition on the lens parameters for this metamorphosis is the 
analogue in the unsheared case of l !23t . There is always a saddle 
on the x-axis, no matter how large 7 is. However, the merger that 
occurs at the saddle on the x-axis may not be the first such merger 
between the two tangential critical curves if 7 is large, as illustrated 
in the panels of Fig.[20] Eqn l |46t gives the condition for the bottom- 
left panel changing to the bottom-right. Once again, this criterion 
holds approximately if r ci , r C2 are small but nonzero. 



6 CONCLUSIONS 

This paper has provided a detailed investigation into the lensing 
properties of binary galaxies, modelled as two isothermal spheres 
with different velocity dispersions. Such models are useful for de- 
scribing strong lenses in which two close galaxies bend light rays 
emanating from a more distant source. With the addition of external 
shear, the models can also represent groups or clusters dominated 
by a pair of massive galaxies. This is of course a very common 



circumstance, with about a quarter of all strong lenses showing ev- 
idence for lensing by pairs or groups of galaxies. 

O ur paper extends and gene ralizes the previous binary lens 
work o flSchneider & Weissl (fl986), who used two point masses. Bi- 
nary isothermal spheres differ qualitatively from the better-known 
case of two point masses by having additional critical curves 
and caustics (or pseudocaustics in the case of singular isothermal 
spheres), and therefore also different multiple-imaging properties. 
If the isothermal spheres have cores, then the number of images 
can be 1, 3, 5 or 7, depending on the location of the source, the 
separation of the binary and the velocity dispersions. If the isother- 
mal spheres are coreless, then the number of images may be 1,2, 
3, 4 or 5. This lens model is sufficiently simple that many of the 
properties can be derived analytically, or nearly so. We have de- 
scribed the shapes of the caustics and critical curves for the case 
of equal and unequal masses, have studied the principal metamor- 
phoses and have classified the critical isochrone topologies. Finally, 
we presented a brief survey of the effects of the inclusion of exter- 
nal shear and internal flattening. 

For comparison with the real world, it is helpful for us to show 
how our models relate to physical systems. The typical velocity dis- 
persion of a large elliptical galaxy is a ~ 250-300 km s _1 , while 
typical core radii are in the range ~ 10-100 pc. Given a typical 
lens redshift of z\ ~ 0.5 and source redshift of z s ~ 2, then the 
corresponding Einstein radius of an isolated isothermal sphere is 
~ 10 kpc. The effects studied in this paper become important if 
two galaxies are separated by < 4 Einstein radii (see e.g., Fig [TJ. 
At the typical redshifts, this corresponds to a separation of < 40 
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Figure 19. The lens here is two cored isothermal spheres at (±a, 0) with external shear of 7 = 0.2 in the direction <f> = 5°. Without shear the critical 
curves and caustics would be Type 1 . The positions of two point sources are shown with the caustics in the left panel; the right panels are contours of the 
Fermat potential for the two sources, together with saddle-point images and their critical isochrones (coloured points and lines) and the images (black points) 
corresponding to Fermat maxima and minima. The order of arrival of the saddle-point images is red, yellow, green, blue, magenta. The four saddles of the 
nine-image case give a critical isochrone topology Li+ Li+ Le~ Le~, and two of the images are probably unobservable 'core' images; the five saddles of 
eleven-image case give Li + Li+ Le~ Le~ Le~ and nine of the images are observable. 



kpc (or < 5"). A good example is provided by CASSOWARY 5, 
which has two lensing galaxies separated by ~ 5" l lBelokurov et alj 
2008). The three images straddle the two lensing galaxies, much as 
the illustration in the lower panel of Fig.QJJ] CASSOWARY 5 looks 
like an example of our image configuration 3B-1, though it impor- 
tant to bear in mind that the existing, primarily Sloan Digital Sky 
Survey, data is comparatively shallow (r < 23). Deeper follow-up 
imaging may revel the presence of more images in CASSOWARY 
5. 

There are some interesting ways to extend the work in this pa- 
per. First, the effects of nearby galaxies are often modelled with 
external shear, so a thorough comparison of the properties of bi- 
nary isothermal lenses with the single isothermal lens plus shear 
is warranted. As is well-known, the Chang-Refsdal lens (or point 
mass plus external shear) is only partially successful in describ- 
ing the lensing properties of binary point masses. Second, although 
we have sketched the principal effects of flattening and shear here, 
it would be interesting to prove the maximum number of possible 
images. The fact that caustics can have points of self-intersection 
means that regions of the source plane are created for which very 
high-order imaging is possible. Third, there are still very few highly 
non-axisymmetric lenses that have been studied. The binary lens 
(whether point masses or double galaxies) is one example, but it 
would be interesting to study the lensing properties of needles, fil- 
aments, and sheets, as well. These may have applications to highly 
prolate or highly flattened galaxies, which are the natural end- 
poi nts of the Newto nian collapse of spheroids of collisionless mat- 
ter dLinetalJll965l) . 
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Figure 20. Critical curves of lenses with two singular isothermal spheres 
(dots at (±a, 0)) in very strong external shear. In all four panels, the ex- 
ternal shear is 7 = 0.5, <f> = 10°. The strong shear distorts the tangential 
caustics of the singular isothermal spheres. In particular, it prevents them 
merging on the x-axis until Ei are much larger: without external shear, 
the lenses in all four panels would have a single common tangential caus- 
tic. The condition {46) holds at the transition between the topologies of the 
bottom-left and bottom-right panels. 
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Figure 21. The critical curves (left panel) and caustics (right panel) of a 
double singular isothermal sphere lens with a = 10, E\ = 3.0, E% = 6.0 
and moderate external shear (7 = 0.2) in the direction <f> = 40°, which 
distorts the tangential critical curves (and caustics). It is believable, but not 
obvious, that the critical curves will meet on the x-axis if Ei are increased. 
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